Intermittency corrections to the mean square particle acceleration in high Reynolds 

number turbulence 
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The mean square particle acceleration in high Reynolds number turbulence is dominated by 
the mean square pressure gradient. Recent experiments by Voth et al [Phys. Fluids 10, 2268 
(1998)] indicate that this quantity, when normalized by the 1941 Kolmogorov value, e 3//2 ^~ 1//2 , is 
independent of Reynolds number at high Reynolds number. This is to be contrasted with direct 
numerical simulations of Vedula and Yeung [Phys. Fluids 11, 1208, (1999)] which show a strong 
increase of the same quantity with increasing Reynolds number. In this paper we suggest that there 
is no inherent conflict between these two results. A large part of the increase seen in DNS is shown 
to be associated with finite Reynolds number corrections within the quasi-Gaussian approximation. 
The remaining intermittency corrections are relatively slowly increasing with Reynolds number, and 
are very sensitive to subtle cancellations among the longitudinal and transverse contributions to the 
mean square pressure gradient. Other possible theoretical subtleties are also briefly discussed. 
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There has been considerable recent interest in the 
study of fluid particle acceleration in turbulent flows. A 
natural focus for study is the scaled mean square fluid 
particle acceleration 

/? = i/VV^V}, (1) 

where a = dv/dt is the fluid particle acceleration, v is 
the kinematic viscosity of the fluid, and e is the average 
rate of energy dissipation per unit mass. In the 1941 Kol- 
mogorov theory, (3 is independent of Reynolds number. 

Direct numerical simulations of isotropic turbulence by 
Vedula and Yeung jlj and by Gotoh and Rogallo ||] show 
that the scaled acceleration (|j) increases strongly with 

1/2 

Reynolds number, approximately as i? A , where R\ is 
the Taylor microscale Reynolds number. These DNS also 
give the expected result that the dominant contribution 
to ft comes from the mean square pressure gradient. The 
viscous contribution is no more than a few percent, and 
will not be considered further in this paper. The DNS are 
currently limited to a maximum R\ of 235. By contrast, 
a recent high resolution measurement of fluid particle 
acceleration by Voth et al || suggests that j3 is nearly 
independent of Reynolds number in a range of R\ be- 
tween 1000 and 2000. This result is for the flow between 
counter-rotating disks. The quantity /3 used in this pa- 
per is a factor of 3 times the scaled acceleration variance 
ao used in [|| and Q. These references look at a sin- 
gle Cartesian component, and we look at the sum over 
Cartesian components. 

We start from an exact result for the mean square 
pressure gradient for isotropic turbulence due to Hill and 
Wilczak M, 

/>oo 

((Vp) 2 )=4/ r- 3 [L(r)+T(r)-6M(r)]dr, (2) 
Jo 



where we have taken the density p = 1 with no loss of 
generality. In (||) the quantities L{r) 7 T(r), and M(r) 
are the three independent fourth order velocity structure 
functions. We refer to these as the longitudinal L(r), 
transverse T(r), and mixed M(r), and they are given by 

L{r) = (A« 4 (r)>, T(r) = (A« 4 (r)), M(r) = (Au 2 {r)Av 2 {r)). 

(3) 

In a recent letter ||, we have analyzed a similar ex- 
pression for the pressure structure function D p (r) in the 
inertial range, and we found it to be very sensitive to 
strong cancellations among the positive terms propor- 
tional to L(r) and T(r) and the negative term propor- 
tional to M{r). We expect a comparable degree of can- 
cellation here for the mean square pressure gradient. To 
estimate this cancellation, it is useful to introduce the 
quasi-Gaussian approximation, where 

L{r) = i[D LL {r)] 2 ,T(r) = 3[D N N {r)} 2 , M (r) = D LL (r)D NN ( 

(4) 

In (||), Dll{t) is the second order longitudinal veloc- 
ity structure function, and D n n (r) is the second order 
transverse velocity structure function. For isotropic tur- 
bulence, these two quantities are related by 

D NN (r) = D LL {r) + (r/2)[dD LL (r)/dr}. (5) 

If (|j) and (|5|) are substituted into we obtain 

/>oo 

<(Vp) 2 )=3/ r-^dDLL^/drfdr, (6) 
Jo 
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which was derived by Obukhov in 1949 || and by Batch- 
elor in 1951 jjj. To evaluate (^) at high Reynolds num- 
bers, Batchelor introduced an interpolation formula for 
the longitudinal structure function, which remains an ac- 
curate and useful approximation today. He suggested 
that 

Dll{t) = (l/15)e 2 / V 2 / 3 x 2 (l + ax 2 y 2 '\ (7) 
where 

x = r/r), 

e is the average rate of energy dissipation per unit mass, 
r\ = j/ 3 / 4 e~ 1 / 4 , is the Kolmogorov dissipation length 
scale, and v is the kinematic viscosity of the fluid. The 
adjustable parameter a is taken asa = 0.006455, as dis- 
cussed elsewhere [||. Substituting (0) into (||), the inte- 
gral in (H) is easily evaluated to give 

pQG = £-3/2^1/2^)2^ = ( 4 / 175a ) w 3 5 ( 8 ) 

A survey of earlier results for the value of (3® G is given 
by Kaneda || . 

At low Reynolds number, there are two distinct correc- 
tions to (||) to be expected. Within the quasi-Gaussian 
approximation, there will be a finite Reynolds number 
correction due to the cutoff of (||) at large r. This occurs 
approximately at the integral length scale L. This leads 
to a viscosity independent subtraction from the mean 
square pressure gradient in (^]). Since the high Reynolds 
number limit of ((Vp) 2 } is proportional to v" 1 ' 2 ~ R\, 
the quasi-Gaussian value of j3 should have a Reynolds 
number dependence, 

^ G - B(l - C/R X ) (9) 

It is difficult to estimate the coefficient C from theory, 
but fortunately the DNS can separately evaluate the full 
value of f3 and its approximate value from the quasi- 
Gaussian approximation. These values, along with the 
ratio of the full to quasi-Gaussian contribution are given 
in Table 1. The full value is from [Q, and the quasi- 
Gaussian value has been computed from the same data 
by Vedula and Yeung, and sent to us privately. The full 
value is three times the value reported in table 1 of 
since it is the sum over all three Cartesian components 
which is used here. 

The values of (3 QG are well fit by (§) with B = 3.75 
and C = 14.7. The value of B is reasonably close to the 
value of B = 3.5 in (|J). We have no interpretation for 
the numerical value of C . Since the large r cutoff is in 
the energy containing range, this value should depend on 
the large scale flow conditions, and is not expected to be 
universal. 

The ratio of the full value to the quasi-Gaussian value 
is determined by non-Gaussian corrections in the dissipa- 
tion range of scales r. These corrections are well known 
to increase as the scale size r decreases. In the incrtial 



range, they are characterized by the well-known anoma- 
lous scaling exponents for the fourth order structure func- 
tions. The situation in the dissipation range is more com- 
plicated and less well characterized, but the dominant 
physical phenomenon is still the increasing intermittency 
with decreasing scale size. We thus refer to the correc- 
tions to the quasi-Gaussian value as intermittency correc- 
tions. In Table 1, we see that /3^ G increases by a factor 
of 1.58 from the smallest to the largest Reynolds num- 
bers in the DNS, and that the intermittency correction 
increases by a factor 1.66 in the same range. For higher 
Reynolds numbers, however, (^) indicates that fi^ G is es- 
sentially constant. The Reynolds number dependence of 
the intermittency correction is difficult to estimate since 
the cancellation of the terms in (||) can be a sensitive 
function of Reynolds number. From Table 1, these cor- 
rections appear to increase approximately as 

f3/p^ G = Const.R Q x 23 , (10) 

but the power law fit is not very accurate, showing 
marginal signs of flattening at higher R\. 

The experiment by Voth et al || is in the flow be- 
tween counter-rotating disks in a range of R\ from 1000 
to 2000. The measured value of (3 in this range is inde- 
pendent of R\ and is 21 ± 9. Most of the uncertainty 
is in the absolute value. The Reynolds number depen- 
dence is known fairly accurately. If we extrapolate the 
intermittency correction from the DNS of Vedula and 
Yeung, we obtain a value of (3 — 14.0 at R\ = 1000 and 
(3 = 16.4 at R x = 2000. We do not know if the exper- 
iment is accurate enough to exclude this weak increase 
in this Reynolds number range, but the absolute value 
is large enough that substantial corrections to the quasi- 
Gaussian approximation are surely present. 

Now let us examine some of the subtle effects which can 
cause a complicated Reynolds number dependence of the 
intermittency corrections. First look at (|^) in the quasi- 
Gaussian approximation. All of the integrands show a 
broad peak in the neighborhood of r = IO77. Substitut- 
ing (H), (||) and (||) into (||), all of the integrals can be 
done analytically to give 

[3 QG = (12/225a)(1.5000 + 3.4286 - 4.5000) = (4/175a) « 3.5, 

(11) 

where the first positive term comes from L(r), the sec- 
ond positive term from T(r), and the negative term from 
M(r). The cancellation among the three terms is sub- 
stantial. If the intermittency corrections to each term 
have slightly different Reynolds number dependence, this 
slight difference can be magnified in an unpredictable 
way in the net value of (3. Further, the Hill-Wilczak 
formula (||) is for isotropic turbulence, and is sensitive 
to corrections due to anisotropy. The experiment is not 
fully isotropic even at dissipation scales. Corrections due 
to anisotropy can be substantial, and nothing is known 
about their Reynolds number dependence. 
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Neglecting such possible differences, it is still a subtle 
and difficult question to estimate the Reynolds number 
dependence of the intermittency correction to the lon- 
gitudinal contribution. This is governed, to a good ap- 
proximation by F(10r))/3, where F(r) is the longitudinal 
flatness factor 

F(r) = <A U 4 (r))/(A W 2 (r)) 2 , (12) 

which is 3 for a Gaussian distribution. The crossover of 
this flatness from the inertial to the dissipation range is 
quite complicated because of the multifractal distribu- 
tion of dissipation length scales. A specific model of this 
multifractal crossover has been calculated by Eggers and 
Wang jicj]. Working with the input data to Fig. 6 of that 
paper, we have calculated the correction factor F(l0rj)/3 
versus Reynolds number in the Reynolds number range 
of H . It increases somewhat more strongly than the ex- 
trapolation from the DNS, and thus does not by itself 
explain the observed plateau in the scaled acceleration 
variance. 

We note that the longitudinal velocity derivative flat- 
ness in the flow between counter-rotating disks has been 
observed jll[] to have a plateau in its Reynolds number 
dependence in the same range of Reynolds numbers stud- 
ied by Voth et al. The origin of this plateau is not un- 
derstood. It is not known if it is a general feature of 
turbulence at these Reynolds numbers, or if it is a spe- 
cial feature of the flow between counter-rotating disks. 
The present study started with the idea that these two 
experiments might be related, but we have not found any 
simple connection between them. 

To summarize, there is no essential disagreement be- 
tween the low Reynolds number DNS and the high 
Reynolds number experiment. In the DNS, a large part 
of the increase in the scaled mean square pressure gra- 
dient (5 arises from finite Reynolds number corrections 
within the quasi-Gaussian approximation, and these cor- 
rections are obtained directly from the DNS. These cor- 
rections are negligible for the high Reynolds number ex- 
periment. The remaining intermittency corrections are 
increasing rather slowly with Reynolds number, and it is 
very difficult to estimate how they will extrapolate to the 
Reynolds number range of the experiment. The experi- 
mental value of (3 — 21 ± 9 has considerable uncertainty 
in its absolute value, but it is almost surely much larger 
than the quasi-Gaussian value, which can be accurately 
calculated as (3® G ~ 3.5. The approximate Reynolds 
number independence observed for j3 does not mean that 
the experiment is consistent with 1941 Kolmogorov scal- 
ing, but only that the intermittency corrections, though 
large, are not varying rapidly with Reynolds number in 
the range of the experiments. 
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3.60 
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1.60 
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6.81 


3.19 


2.14 


140 


8.16 


3.37 


2.42 


235 


9.39 


3.55 


2.65 



TABLE I. Scaled mean square pressure gradient and its 
quasi-Gaussian value from the DNS of Vedula and Yeung. 
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